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Rescaled Perturbation Theory 

Tomoya Hayata*) 

Department of Physics, The University of Tokyo, Tokyo 113-0033, Japan 

A nonperturbative method that can go beyond the weak coupling perturbation theory 
is introduced. The essential idea is to formulate a set of exact differential equations as a 
function of coupling strength g. Unlike other resummation methods in which information 
on higher-order terms is necessary, we only need a leading-order perturbative formula in 
every step to reach a large value of g. The method is applied to the quantum anharmonic 
oscillator and quantum double-well potential in one dimension. Both are known to have 
divergent series in the weak coupling perturbation and the latter is not Borel-summable. 
Our method is shown to work well from the weak coupling to the strong coupling for the 
energy eigenvalues and the wave functions. The method is also applied successfully to a 
D ■ system with a time-dependent external field. 

In quantum many-body problems and in quantum field theories, there are many 
■ systems where one needs to study phenomena over a broad range of coupling strengths.^ 

Typical examples are the BEC-BCS crossover in ultracold atoms where the coupling 
strength is changed by utilizing the Feschbach resonance hadron-quark transi- 
tions in dense QCD (quantum choromodynamics) where the coupling runs with the 
baryon density® and the quark-gluon plasma near the critical temperature where 
the coupling runs with temperature.-^ Thus far, a number of techniques have been 
developed to resum naive perturbative series or to produce a better convergent se- 
ries under the names of resummed perturbation, optimized perturbation, variational 
perturbation and so on. (For reviews, see e.g., Refs. ED and [7]).) 
CN ■ In this article, we introduce a novel method that can go beyond the weak coupling 

perturbation theory: we test its applicability to quantum mechanical examples such 
| as the anharmonic oscillator (AHO) and double- well potential (DWP). Instead of 

resumming the perturbation series, we iterate the lowest order perturbation through 
the rescaling of the coupling strength, so that the method gives a global coupling 
dependence of eigenvalues and eigenvectors simultaneously. The basic formula has 
formal resemblance to the exact renormalization group method (ERG, see e.g., Ref. 
8) for a review) if we replace the cutoff scale of ERG by the coupling strength in the 
present approach. 

Let us consider a quantum system such that the Hamiltonian is 

£ : n( g ) = n + gn int , (i) 

where g is a dimensionless coupling strength and Ho is assumed to be solved accu- 
rately. Our goal is to obtain E n (g) and \ijj n (g)} satisfying 

n( 9 )\Mg)) = E n (g)\M9)) ■ (2) 

To simplify the discussion, we consider a nondegenerate %q. Generalization to the 
degenerate case is straightforward. Suppose we know E n (g) and \ip n (g)}, then we can 



*' E-mail: t.hayata@nt.phys.s.u-tokyo. ac.jp 



2 



Letters 



calculate E n {g + 8g) and \ip n {g + Sg)) as long as 5g <C 1 by using the leading-order 
Rayleigh-Schrodinger perturbation theory (RSPT) for T~L(g + <5g) =1-L(g) + SgHmt, 



E n (g + Sg) = E n (g) + 5g{ip n (g)\n int \Mg)) , (3) 

^(g)l^intl^n(g] 
£„(<?) - Ei( 9 ) 



\Mg + Sg)) = \Mg)) + s 9 Y J ^f^if Mg)) • (4) 



Then, we can set up differential equations for the eigenvalues, matrix elements 
H-^ t (g) = (tpi(g)\'Hint\tpj(g)} and the wave functions as 

^Ekig) = Bgt(g) , (5) 

±tfi (n) = V g£Cg)g£(g) + V H int(g) H ^(g) rfil 
^ intW Ei {g) - E k (g) + j£ E {g) - E k (g) ' { > 

Note that we have no approximation to obtain the right-hand sides of these equations, 
so that they are exact equations. We take {ipi(0)\ipi(g)) = (ipi(0)\ipi(0)) = 1 as a 
normalization of the state vector, so that the norm of eigenvectors is conserved, i.e., 
^(i ; i(g)\' l l ; i(g)) = 0- Coupled differential equations ([5]), © and ([7]) can be solved 
under the initial conditions, Ei(0), HjjL(O) and \ipi(Q)) and one may go to arbitrary 
large values of g in principle. 

Let us now apply this method to the one-dimensional quantum AHO: 

^ahoG?) = \p 2 + \x 2 + gx 4 . (8) 

This is a typical example that RSPT leads to a divergent series but is Borel- 
summable. (See e.g., Refs. [6]) and [9|) for reviews.) Using Cij(g) = (ipj(0)\'ipi(g)) , 
we can rewrite Eq. ([7]) as 

The initial conditions at g = read, 

Ei(0)=i + l/2 , (10) 

H£ t (0) = 1 [ dx e -^x*fli(x)^(x) , (11) 

c ij (0) = 5 ij , (12) 

where Hi{x) is the i-th Hermite polynomial. In actual calculation, we solve the 
coupled differential equations for the states i = 0, 1, • • • N — 1. For g up to 10, 
taking N = 50 is accurate enough to have the eigenvalues of the low-lying states as 
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Table I. Comparison of the lowest three eigenvalues of AHO. The upper numbers are obtained from 
the present method and the lower numbers from Ref. I10[) . 
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0.1 


0.55914633 


1.7695026 


3.1386243 




0.55914633 


1.7695026 


3.1386243 


0.5 


0.69617582 


2.3244064 


4.3275250 




0.69617582 


2.3244064 


4.3275250 


1.0 


0.80377065 


2.7378923 


5.1792917 




0.80377065 


2.7378923 


5.1792915 


5.0 


1.2245874 


4.2995081 


8.3179758 




1.2245870 


4.2995017 


8.3179605 
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1.5049814 


5.3216308 


10.348359 




1.5049724 


5.3216080 


10.347056 




shown in Table HI Here, we compare our results with those given in Ref. [10|) by the 
diagonalization method. 

Next, we study the quantum DWP with the Hamiltonian, 

n DWP (g) = l -v 2 - gx 2 + \x 4 , (g > 0) (13) 

where the RSPT also breaks down and is known to be Borel-nonsummable. (See 
e.g., Refs. [6]) and [9]) for reviews.) To make a firm connection between DWP and 
AHO in our approach, we rewrite the Hamiltonian as 

Ui DWP (9')=(\p 2 + \x 2 + \x i )-g'x 2 (14) 

= %AHO - Spirit , (15) 

so that the unperturbed Hamiltonian at g' = becomes AHO. Note that the system 
changes from AHO for g' < 1/2 to DWP for g' > 1/2. In Fig. [IJ the shape of the 
potential is shown for several different values of g' . 
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Table II. Comparison of the lowest two eigenvalues for the DWP. The upper numbers are calculated 
by our method, while the lower numbers are taken from Ref. 
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0.50 


0.53018104538 


1.8998365150 




0.53018104524 


1.8998365149 


1.0 


0.32882650295 


1.4172681012 




0.32882650260 


1.4172681011 


5.5 


-10.316788242 


-10.316773352 




-10.316788351 


-10.316773442 


8.0 


-25.420689499 


-25.420692377 




-25.420693642 


-25.420693642 



Using the results of the previous AHO for N = 50, we calculate the energy 
eigenvalues of this system as a function of g' . The differential equations are the 
same as Eqs. ([5|), ([6|) and ([9]), where iV'j(O)) is the eigenvector of AHO in the present 
case. The initial conditions read, 

EM = Er°(l) , (16) 



,2, 

C t (0) = (^(0)|-x 2 |^(0)) (17) 

N-l 

= E4-( )4r^ 7 J m /*.-■<-.>)*.(.)*.(.> . 

k,i=o v /IZ J 

c lJ (0)=5 ij . (18) 

Here, i,j runs from through N — 1, and N is taken to be 50 as before. This is 
good enough again for low- lying eigenvalues as shown in Table ITU where comparison 
to the high- accuracy eigenvalues given in Ref. 11) by the diagonalization method is 
shown. 

In Fig. [21 we show the g' dependence of the lowest six energy eigenvalues from 
the weak coupling to the strong coupling. As g' increases, the degeneracy of the 
even-n and odd-n eigenstates takes place starting from the low-energy state. This 
can be clearly seen by looking at the probability distributions of the lowest two states 
as shown in Fig. [3J Depending on g' , the energy eigenstate changes from a single 
oscillator to the superposition of double (left and right) oscillators. 

The present method can be generalized to the time-dependent Schrodinger equa- 
tion with a nonadiabatic time-dependent potential: 

id t \m) = nt)\m) , (19) 

where H(t) = 7~Lq + g{t)Tl\ n t. We expand the state vectors by the adiabatic basis: 

\m) = Y. a n(t)Ht)), (20) 

n 

H(t)\n(t)) = E n (t)\n(t)) , (21) 
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Fig. 2. Lowest six energy eigenvalues for the DWP with N = 50. Filled circles are high-accuracy 
numerical results given in Ref. Illj) . 




Fig. 3. Probability distributions for the lowest two states for the DWP with N = 50. The solid 
line is for the ground state and the dashed line is for the first excited state. The lighter gray 
line is for g = 0.5, the gray line is for g' = 2.5 and the black line is for g' = 6.0. 



where E n (t) and \n(t)) are instantaneous eigenvalues and eigenvectors with their 
time dependence implicitly through the coupling g(t). We perform a unitary trans- 
formation: 

a n {t) = a n (t)e~ i0 ^ , (22) 

Q n {t) = [ dt' E n (t') , (23) 
J o 

and obtain 

d t a n (t) = -J2 a m (t)(n(t)\d t \m(t))e i( - e ^- e ^ . (24) 

m 

Although g(t) is an arbitrary smooth function of t, one may always introduce a set of 
time intervals, so that g(t) is a monotonic function of t. Then, within each interval, 
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% = 0.1,_n_=%_ t 



Fig. 4. Lowest three even transition probabilities |a n (t)| 2 (n = 0,2,4) in DWP for N = 50. The 
lighter gray line is for n — 0, the gray line is for n — 2 and the black line is for n — 4. The 
dashed line is for v = 0.1 and the solid line is for v = 3. 



we have 



d g a n {g) - - ^ a m (0)(n(g)|d s |m(g))e 



i(e„(s)-e m (s)) 



■m 

3 1 



d# -E n {g) 

o 5 



(25) 
(26) 



where 5 implies the time derivative of coupling strength. Assuming that %o is n °t 
degenerate, we obtain the following equations from ©, © and ((TJ): 



9 n[9) % E n (g) - E m (gf m[9)e 
d g On(g) = -^E n {g) . 



(27) 
(28) 



To calculate the evolution of a n (g) and n (g), we need the coupling dependence 
of E n (g) and H^™(g), which is directly obtained from our approach, so that we 
can trace the flows of E n (g), a n (g) and n (g) simultaneously through the flows of 

Here, we solve the differential equations numerically in the case that the coupling 
strength depends linearly on time (i.e., g = vt) and obtain the transition probabilities 
to adiabatic states |a n (i)| 2 . In Fig. [H taking the ground state of AHO as an initial 
condition, we show the transition probabilities as a function of time for the DWP with 
the Hamiltonian: 'H^plff'fi) = vt) in Eq. ((HJ. Because of the parity conservation, 
the transition probability for odd n becomes zero. If we change the potential slowly 
as the dashed lines (y = 0.1) in Fig. [H the instantaneous basis at n = is mostly 
occupied at any time. On the other hand, if we change the potential suddenly as the 
solid lines (v = 3.0), higher instantaneous states are excited with mutual oscillations. 

In summary, we introduced a nonperturbative method applicable to large values 
of the coupling strength g. The essential idea is to formulate a set of exact differential 
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equations as a function of g, so that they can be solved under appropriate known 
initial conditions to reach large values of g. Unlike resummation methods in which 
information on higher-order terms in the naive perturbation series is necessary, we 
only need a leading-order perturbative formula in every step. This is similar to the 
situation in the exact renormalization group method. In the present method, all 
the eigenvalues and eigenstates can be calculated simultaneously and accurately as 
a function of g as long as we prepare a sufficiently large enough dimension N of the 
Hilbert space. Note here that we do not need to diagonalize the N x N Hamiltonian; 
instead, we trace the flows of eigenvalues and eigenvectors starting from the solvable 
N x N Hamiltonian. We have applied our method to the quantum anharmonic 
oscillator and quantum double-well potential in one dimension. Both are known to 
have divergent series in RSPT and the latter is not even Borel-summable. We found 
that our method works well from the weak coupling to the strong coupling for the 
energy eigenvalues and wave functions. Furthermore, because of the flow equations 
of the energy eigenvalues and the matrix elements H-^ t (g), we could solve even the 
time-dependent Schrodinger equation for the potential with nonadiabatic variation 
in time. 

The basic idea of the present method can also be generalized to the quantum 
field theory and quantum statistical mechanics. In these cases, one can derive a 
differential equation for the partition function Z[J,g] as functions of the external 
field J and the coupling strength g, 



where the initial condition Z[J, 0] is assumed to be solvable. Similarly, we can also 
apply the idea to the generalization of the linear response theory to derive a formula: 



Here, A is an external field, j is a conjugate quantity and h denotes a Heisenberg 
picture in terms of H[t;g] = %Q+gjA(t). Applications of these equations to physical 
systems will be discussed in forthcoming publications. 

After the completion of this work, we became aware that differential equations 
similar to ours are discussed in Appendix C of Ref. 12), although the equations 
are applied in a rather different context from ours. The author would like to thank 
Tetsuo Hatsuda, Kyogo Kawaguchi and Haruki Watanabe for useful discussions and 
suggestions. 
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